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Abstract 

Measurement-based quantum computation has emerged from the physics community as a 
new approach to quantum computation where the notion of measurement is the main driving 
force of computation. This is in contrast with the more traditional circuit model which is based 
on unitary operations. Among measurement-based quantum computation methods, the recently 
introduced one-way quantum computer |RB01| stands out as fundamental. 

We develop a rigorous mathematical model underlying the one-way quantum computer and 
present a concrete syntax and operational semantics for programs, which we call patterns, and an 
algebra of these patterns derived from a denotational semantics. More importantly, we present 
a calculus for reasoning locally and compositionally about these patterns. We present a rewrite 
theory and prove a general standardization theorem which allows all patterns to be put in a 
semantically equivalent standard form. Standardization has far-reaching consequences: a new 
physical architecture based on performing all the entanglement in the beginning, parallclization 
by exposing the dependency structure of measurements and expressiveness theorems. 

Furthermore we formalize several other measurement-based models e.g. Teleportation, Phase 
and Pauli models and present compositional embeddings of them into and from the one-way 
model. This allows us to transfer all the theory we develop for the one-way model to these mod- 
els. This shows that the framework we have developed has a general impact on measurement- 
based computation and is not just particular to the one-way quantum computer. 

1 Introduction 

The emergence of quantum computation has changed our perspective on many fundamental aspects 
of computing: the nature of information and how it flows, new algorithmic design strategies and 
complexity classes and the very structure of computational models [NCOOj . New challenges have 
been raised in the physical implementation of quantum computers. This paper is a contribution to 
a nascent discipline: quantum programming languages. 

This is more than a search for convenient notation, it is an investigation into the structure, 
scope and limits of quantum computation. The main issues are questions about how quantum 
processes are defined, how quantum algorithms compose, how quantum resources are used and how 
classical and quantum information interact. 

Quantum computation emerged in the early 1980s with Feynman's observations about the dif- 
ficulty of simulating quantum systems on a classical computer. This hinted at the possibility of 
turning around the issue and exploiting the power of quantum systems to perform computational 
tasks more efficiently than was classically possible. In the mid 1980s Deutsch [Deu87] and later 
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Deutsch and Jozsa |DJ92] showed how to use superposition - the ability to produce linear combi- 
nations of quantum states - to obtain computational speedup. This led to interest in algorithm 
design and the complexity aspects of quantum computation by computer scientists. The most 
dramatic results were Shor's celebrated polytime factorization algorithm |Sho94j and Grover's sub- 
linear search algorithm [Gro98j . Remarkably one of the problematic aspects of quantum theory, 
the presence of non-local correlation - an example of which is called "entanglement" - turned out 
to be crucial for these algorithmic developments. 

If efficient factorization is indeed possible in practice, then much of cryptography becomes 
insecure as it is based on the difficulty of factorization. However, entanglement makes it possible 
to design unconditionally secure key distribution [BB841 lEke91| . Furthermore, entanglement led 
to the remarkable - but simple - protocol for transferring quantum states using only classical 
communication [BBC + 93] ; this is the famous so-called "teleportation" protocol. There continues to 
be tremendous activity in quantum cryptography, algorithmic design, complexity and information 
theory. Parallel to all this work there has been intense interest from the physics community to 
explore possible implementations, see, for example, |NC00| for a textbook account of some of these 
ideas. 

On the other hand, only recently has there been significant interest in quantum programming 
languages; i.e. the development of formal syntax and semantics and the use of standard machinery 
for reasoning about quantum information processing. The first quantum programming languages 
were variations on imperative probabilistic la nguag es and emphasized logic and program develop- 



ment based on weakest preconditions [SZ00[ lOOl . The first definitive treatment of a quantum 



programming language was the flowchart language of Selinger |Sel04b| . It was based on combining 
classical control, as traditionally seen in flowcharts, with quantum data. It also gave a denotational 
semantics based on completely positive linear maps. The notion of quantum weakest preconditions 
was developed in [DP06j . Later people proposed languages based on quantum control |AG05| . The 
search for a sensible notion of higher- type computation [SV051 lvT04j continues, but is problem- 
atic [Sel04cj . 

A related recent development is the work of Abramsky and Coecke |AC041 ICoe04j where they 
develop a categorical axiomatization of quantum mechanics. This can be used to verify the correct- 
ness of quantum communication protocols. It is very interesting from a foundational point of view 
and allows one to explore exactly what mathematical ingredients are required to carry out certain 
quantum protocols. This has also led to work on a categorical quantum logic [AD04| . 

The study of quantum communication protocols has led to formalizations based on process 
algebras [GN051 IJL04] and to proposals to use model checking for verifying quantum protocols. A 
survey and a complete list of references on this subject up to 2005 is available |Gay05| . 

These ideas have proven to be of great utility in the world of classical computation. The use of 
logics, type systems, operational semantics, denotational semantics and semantic-based inference 
mechanisms have led to notable advances such as: the use of model checking for verification, 
reasoning compositionally about security protocols, refinement-based programming methodology 
and flow analysis. 

The present paper applies this paradigm to a very recent development: measurement-based 
quantum computation. None of the cited research on quantum programming languages is aimed 
at measurement-based computation. On the other hand, the work in the physics literature does 
not clearly separate the conceptual layers of the subject from implementation issues. A formal 
treatment is necessary to analyze the foundations of measurement-based computation. 
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So far the main framework to explore quantum computation has been the circuit model |Deu89j . 
based on unitary evolution. This is very useful for algorithmic development and complexity analysis 
|BV97] . There are other models such as quantum Turing machines [Deu85j and quantum cellular 
automata |Wat95l lvD96| IDS961 ISW04] , Although they are all proved to be equivalent from the 
point of view of expressive power, there is no agreement on what is the canonical model for exposing 
the key aspects of quantum computation. 

Recently physicists have introduced novel ideas based on the use of measurement and entangle- 
ment to perform computation [GC991 iRB Ol , RBB03, Nic03]. This is very different from the circuit 
model where measurement is done only at the end to extract classical output. In measurement-based 
computation the main operation to manipulate information and control computation is measure- 
ment. This is surprising because measurement creates indeterminacy, yet it is used to express 
deterministic computation defined by a unitary evolution. 

The idea of computing based on measurements emerged from the teleportation protocol [BBC + 93] . 
The goal of this protocol is for an agent to transmit an unknown qubit to a remote agent without 
actually sending the qubit. This protocol works by having the two parties share a maximally en- 
tangled state called a Bell pair. The parties perform local operations - measurements and unitaries 

- and communicate only classical bits. Remarkably, from this classical information the second 
party can reconstruct the unknown quantum state. In fact one can actually use this to com- 
pute via teleportation by choosing an appropriate measurement |GC99] . This is the key idea of 
measurement-based computation. 

It turns out that the above method of computing is actually universal. This was first shown 
by Gottesman and Chuang |GC99| who used two-qubit measurements and given Bell pairs. Later 
Nielsen |Nie03j showed that one could do this with only 4-qubit measurements with no prior Bell 
pairs, however this works only probabilistically. Leung [Leu04j improved this to two qubits, but her 
method also works only probabilistically. Later Perdrix and Jorrand |Per031 IP J04j gave the minimal 
set measurements to perform universal quantum computing - but still in the probabilistic setting 

- and introduced the state-transfer and measurement-based quantum Turing machine. Finally 
the one-way computer was invented by Raussendorf and Briegel (RB01, RBQ2] which used only 
single-qubit measurements with a particular multi-party entangled state, the cluster state. 

More precisely, a computation consists of a phase in which a collection of qubits are set up in a 
standard entangled state. Then measurements are applied to individual qubits and the outcomes of 
the measurements may be used to determine further measurements. Finally - again depending on 
measurement outcomes - local unitary operators, called corrections, are applied to some qubits; this 
allows the elimination of the indeterminacy introduced by measurements. The phrase "one-way" 
is used to emphasize that the computation is driven by irreversible measurements. 

There are at least two reasons to take measurement-based models seriously: one conceptual 
and one pragmatic. The main pragmatic reason is that the one-way model is believed by physicists 
to lend itself to easier implementations [Nie04l ICAJ051 IBR051 ITPKV041 ITPKV061 IWkJRR+05l 
IKPA061 IBES051 ICCWD061 IBBFM06] . Physicists have investigated various properties of the cluster 
state and have accrued evidence that the physical implementation is scalable and robust against 
decoherence |Sch03l IHEB041 IDAB031 ldNDM04bl ldNDM04al IMP041 IGHW051 IHDB051 IDHN06] . 
Conceptually the measurement-based model highlights the role of entanglement and separates the 
quantum and classical aspects of computation; thus it clarifies, in particular, the interplay between 
classical control and the quantum evolution process. 

Our approach to understanding the structural features of measurement-based computation is to 
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develop a formal calculus. One can think of this as an "assembly language" for measurement-based 
computation. Ours is the first programming framework specifically based on the one-way model. We 
first develop a notation for such classically correlated sequences of entanglements, measurements, 
and local corrections. Computations are organized in patterns^, and we give a careful treatment 
of the composition and tensor product (parallel composition) of patterns. We show next that such 
pattern combinations reflect the corresponding combinations of unitary operators. An easy proof 
of universality follows. 

So far, this is primarily a clarification of what was already known from the series of papers 
introducing and investigating the properties of the one-way model |RB01l |RB02| |RBB03| . However, 
we work here with an extended notion of pattern, where inputs and outputs may overlap in any 
way one wants them to, and this results in more efficient - in the sense of using fewer qubits - 
implementations of unitaries. Specifically, our universal set consists of patterns using only 2 qubits. 
From it we obtain a 3 qubit realization of the R z rotations and a 14 qubit realization for the 
controlled-?/ family: a significant reduction over the hitherto known implementations. 

The main point of this paper is to introduce a calculus of local equations over patterns that 
exploits some special algebraic properties of the entanglement, measurement and correction op- 
erators. More precisely, we use the fact that that 1-qubit XY measurements are closed under 
conjugation by Pauli operators and the entanglement command belongs to the normalizer of the 
Pauli group; these terms are explained in the appendix. We show that this calculus is sound in that 
it preserves the interpretation of patterns. Most importantly, we derive from it a simple algorithm 
by which any general pattern can be put into a standard form where entanglement is done first, 
then measurements, then corrections. We call this standardization. 

The consequences of the existence of such a procedure are far-reaching. Since entangling comes 
first, one can prepare the entire entangled state needed during the computation right at the start: 
one never has to do "on the fly" entanglements. Furthermore, the rewriting of a pattern to stan- 
dard form reveals parallelism in the pattern computation. In a general pattern, one is forced to 
compute sequentially and to strictly obey the command sequence, whereas, after standardization, 
the dependency structure is relaxed, resulting in lower computational depth complexity. Last, the 
existence of a standard form for any pattern also has interesting corollaries beyond implementation 
and complexity matters, as it follows from it that patterns using no dependencies, or using only the 
restricted class of Pauli measurements, can only realize a unitary belonging to the Clifford group, 
and hence can be efficiently simulated by a classical computer |Got97] . 

As we have noted before, there are other methods for measurement-based quantum comput- 
ing: the teleportation technique based on two-qubit measurements and the state-transfer approach 
based on single qubit measurements and incomplete two-qubit measurements. We will analyze the 
teleportation model and its relation to the one-way model. We will show how our calculus can be 
smoothly extended to cover this case as well as new models that we introduce in this paper. We 
get several benefits from our treatment. We get a workable syntax for handling the dependencies of 
operators on previous measurement outcomes just by mimicking the one obtained in the one-way 
model. This has never been done before for the teleportation model. Furthermore, we can use this 
embedding to obtain a standardization procedure for the models. Finally these extended calculi 
can be compositionally embedded back in the original one-way model. This clarifies the relation 
between different measurement-based models and shows that the one-way model of Raussendorf 

1 We use the word "pattern" rather than "program", because this corresponds to the commonly used terminology 
in the physics literature. 
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and Briegel is the canonical one. 

This paper develops the one-way model ab initio but certain concepts that the reader may be 
unfamiliar with: qubits, unitaries, measurements, Pauli operators and the Clifford group are in 
an appendix. These are also readily accessible through the very thorough book of Nielsen and 
Chuang [NCOOj . 

In the next section we define the basic model, followed by its operational and denotational 
semantics, for completeness a simple proof of universality is given in section [H this has appeared 
earlier in the physics literature |DKP05j . in section [5] we develop the rewrite theory and prove the 
fundamental standardization theorem. In section [6] we develop several examples that illustrate the 
use of our calculus in designing efficient patterns. In section [7] we prove some theorems about the 
expressive power of the calculus in the absence of adaptive measurements. In section [8] we discuss 
other measurement-based models and their compositional embedding to and from the one-way 
model. In section [9] we discuss further directions and some more related work. In the appendix we 
review basic notions of quantum mechanics and quantum computation. 

2 Measurement Patterns 

We first develop a notation for 1-qubit measurement based computations. The basic commands 
one can use in a pattern are: 

• 1-qubit auxiliary preparation iVj 

• 2-qubit entanglement operators E^ 

• 1-qubit measurements M" 

• and 1-qubit Pauli operators corrections X i and Z i 

The indices i, j represent the qubits on which each of these operations apply, and a is a 
parameter in [0, 2ir]. Expressions involving angles are always evaluated modulo 2ir. These types of 
command will be referred to as N, E, M and C. Sequences of such commands, together with two 
distinguished - possibly overlapping - sets of qubits corresponding to inputs and outputs, will be 
called measurement patterns, or simply patterns. These patterns can be combined by composition 
and tensor product. 

Importantly, corrections and measurements are allowed to depend on previous measurement 
outcomes. We shall prove later that patterns without these classical dependencies can only realize 
unitaries that are in the Clifford group. Thus, dependencies are crucial if one wants to define a 
universal computing model; that is to say, a model where all unitaries over <g) n C 2 can be realized. 
It is also crucial to develop a notation that will handle these dependencies. This is what we do 
now. 

2.1 Commands 

Preparation Ni prepares qubit i in state |+).;. The entanglement commands are defined as Eij := 
AZij (controlled-^), while the correction commands are the Pauli operators X { and Z^ 
Measurement Mf is defined by orthogonal projections on 

|+«> := ^(|0> + e-|l» 
|-«> := d?(|0>-e fa |l» 
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followed by a trace-out operator. The parameter a G [0, 2tt] is called the angle of the mea- 
surement. For a = 0, a = ?, one obtains the X and Y Pauli measurements. Operationally, 
measurements will be understood as destructive measurements, consuming their qubit. The out- 
come of a measurement done at qubit i will be denoted by Sj G 7L<i- Since one only deals here with 
patterns where qubits are measured at most once (see condition (Dl) below), this is unambiguous. 
We take the specific convention that Si = if under the corresponding measurement the state 
collapses to \+ a ), and Si = 1 if to |— a ). 

Outcomes can be summed together resulting in expressions of the form s = Yliei s i which we 
call signals, and where the summation is understood as being done in 7L<i- We define the domain 
of a signal as the set of qubits on which it depends. 

As we have said before, both corrections and measurements may depend on signals. Depen- 
dent corrections will be written Xf and Zf and dependent measurements will be written f [M"] s , 
where s,t £7*2 and a G [0, 2tt]. The meaning of dependencies for corrections is straightforward: 
X® = Zf = /, no correction is applied, while X\ = X i and Z\ = Z,-. In the case of dependent 
measurements, the measurement angle will depend on s, t and a as follows: 

*[Mf] s := M^ 1)Sa+tn (1) 

so that, depending on the parities of s and t, one may have to modify the a to one of —a, a + tt 
and —a + tt. These modifications correspond to conjugations of measurements under X and Z: 

XiMfXi = Mr a (2) 
ZiMfZi = Mf +T (3) 

accordingly, we will refer to them as the X and Z-actions. Note that these two actions commute, 
since —a + tt = —a — tt up to 2tt, and hence the order in which one applies them does not matter. 

As we will see later, relations ([2]) and ([3]) are key to the propagation of dependent corrections, 
and to obtaining patterns in the standard entanglement, measurement and correction form. Since 
the measurements considered here are destructive, the above equations actually simplify to 

M?Xi = M~ a (4) 
MfZi = M?- w (5) 

Another point worth noticing is that the domain of the signals of a dependent command, be it a 
measurement or a correction, represents the set of measurements which one has to do before one 
can determine the actual value of the command. 

We have completed our catalog of basic commands, including dependent ones, and we turn 
now to the definition of measurement patterns. For convenient reference, the language syntax is 
summarized in Figure [TJ 

2.2 Patterns 

Definition 1 Patterns consists of three finite sets V, I, O, together with two injective maps l : 
I — > V and o : O — > V and a finite sequence of commands A n . . . A\, read from right to left, applying 
to qubits in V in that order, i.e. A\ first and A n last, such that: 

(DO) no command depends on an outcome not yet measured; 
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0, 1, Si, S + S 



Signals 
Preparations 
Entanglements 
Measurements 
Corrections 



Figure 1: 1-qubit based measurement language syntax 



(Dl) no command acts on a qubit already measured; 

(D2) no command acts on a qubit not yet prepared, unless it is an input qubit; 
(D3) a qubit i is measured if and only if i is not an output. 

The set V is called the pattern computation space, and we write $)v for the associated quantum 
state space (8)jeyC 2 . To ease notation, we will omit the maps i and o, and write simply /, O instead 
of and o(0). Note, however, that these maps are useful to define classical manipulations of 
the quantum states, such as permutations of the qubits. The sets /, O are called respectively the 
pattern inputs and outputs, and we write S)i, and $)o for the associated quantum state spaces. The 
sequence A n . . . A\ is called the pattern command sequence, while the triple (V, /, O) is called the 
pattern type. 

To run a pattern, one prepares the input qubits in some input state ip £ Sji, while the non-input 
qubits are all set to the |+) state, then the commands are executed in sequence, and finally the result 
of the pattern computation is read back from outputs as some (ft £ ¥)o ■ Clearly, for this procedure 
to succeed, we had to impose the (DO), (Dl), (D2) and (D3) conditions. Indeed if (DO) fails, then 
at some point of the computation, one will want to execute a command which depends on outcomes 
that are not known yet. Likewise, if (Dl) fails, one will try to apply a command on a qubit that 
has been consumed by a measurement (recall that we use destructive measurements). Similarly, if 
(D2) fails, one will try to apply a command on a non-existent qubit. Condition (D3) is there to 
make sure that the final state belongs to the output space fioi *.e., that all non-output qubits, and 
only non-output qubits, will have been consumed by a measurement when the computation ends. 

We write (D) for the conjunction of our definiteness conditions (DO), (Dl), (D2) and (D3). 
Whether a given pattern satisfies (D) or not is statically verifiable on the pattern command se- 
quence. We could have imposed a simple type system to enforce these constraints but, in the 
interests of notational simplicity, we chose not to do so. 

Here is a concrete example: 



with computation space {1,2}, inputs {1}, and outputs {2}. To run 7i, one first prepares the first 
qubit in some input state ip, and the second qubit in state |+), then these are entangled to obtain 
AZi2(v^i |+)2)- Once this is done, the first qubit is measured in the |+), |— ) basis. Finally an X 
correction is applied on the output qubit, if the measurement outcome was s± = 1. We will do this 
calculation in detail later, and prove that this pattern implements the Hadamard operator H. 



n 



({l,2},{l},{2},X 2 Sl M°£ 12 iV 2 ) 
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In general, a given pattern may use auxiliary qubits that are neither input nor output qubits. 
Usually one tries to use as few such qubits as possible, since these contribute to the space complexity 
of the computation. 

A last thing to note is that one does not require inputs and outputs to be disjoint subsets of 
V. This, seemingly innocuous, additional flexibility is actually quite useful to give parsimonious 
implementations of unitaries |DKP05j . While the restriction to disjoint inputs and outputs is 
unnecessary, it has been discussed whether imposing it results in patterns that are easier to realize 
physically. Recent work [HEB041 IBR05^ ICAJ05] however, seems to indicate it is not the case. 

2.3 Pattern combination 

We are interested in how one can combine patterns in order to obtain bigger ones. 

The first way to combine patterns is by composing them. Two patterns V\ and V 2 may be 
composed if V\ Pi V 2 = 0\ = I 2 . Provided that V\ has as many outputs as V 2 has inputs, by 
renaming the pattern qubits, one can always make them composable. 

Definition 2 The composite pattern V 2 V\ is defined as: 
-V:=V 1 UV 2 ,I = h, = 2 , 

— commands are concatenated. 

The other way of combining patterns is to tensor them. Two patterns V\ and V 2 may be tensored 
if Vi fl V 2 = 0. Again one can always meet this condition by renaming qubits in a way that these 
sets are made disjoint. 

Definition 3 The tensor pattern V\ % V 2 is defined as: 
- V = Vi U V 2l I = h U h, and O = O x U 2 , 

— commands are concatenated. 

In contrast to the composition case, all the unions involved here are disjoint. Therefore commands 
from distinct patterns freely commute, since they apply to disjoint qubits, and when we say that 
commands have to be concatenated, this is only for definiteness. It is routine to verify that the 
definiteness conditions (D) are preserved under composition and tensor product. 

Before turning to this matter, we need a clean definition of what it means for a pattern to 
implement or to realize a unitary operator, together with a proof that the way one can combine 
patterns is reflected in their interpretations. This is key to our proof of universality. 

3 The semantics of patterns 

In this section we give a formal operational semantics for the pattern language as a probabilistic 
labeled transition system. We define deterministic patterns and thereafter concentrate on them. 
We show that deterministic patterns compose. We give a denotational semantics of deterministic 
patterns; from the construction it will be clear that these two semantics are equivalent. 

Besides quantum states, which are non-zero vectors in some Hilbert space fjy, one needs a 
classical state recording the outcomes of the successive measurements one does in a pattern. If we 
let V stand for the finite set of qubits that are still active (i.e. not yet measured) and W stands 
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for the set of qubits that have been measured (i.e. they are now just classical bits recording the 
measurement outcomes), it is natural to define the computation state space as: 

S := ^v,w^>v x 7% ■ 

In other words the computation states form a V, P^-indexed family of parrel q, T, where q is a 
quantum state from Sjy an d T is a map from some W to the outcome space Z2. We call this 
classical component F an outcome map, and denote by the empty outcome map in Z^. We will 
treat these states as pairs unless it becomes important to show how V and W are altered during a 
computation, as happens during a measurement. 

3.1 Operational semantics 

We need some preliminary notation. For any signal s and classical state V S ZJ^, such that the 
domain of s is included in W, we take sr to be the value of s given by the outcome map T. That is 
to say, if s = Sj, then sr '■= X)j where the sum is taken in Z2. Also if T S ZJj^, and x £ Z2, 
we define: 

T[x/i](i) = x, F[x/i](j) = r(j) for j / i 

which is a map in l^^ 1 ^ . 

We may now view each of our commands as acting on the state space S, we have suppressed V 
and W in the first 4 commands: 



?,r 




g® |+>i,r 


?,r 




AZijq,T 


g,r 






g,r 






FU{i},W,g,r 


'[M°] s 


v,wu{i},(+ ar \ iq ,r[o/i 


FU{i},W,g,r 


*[Mf] a 


V,WU{i},(- ar \ iq ,T[l/i 



where «r = (— l) Sr a + trft following equation ([T|). Note how the measurement moves an index 
from V to W; a qubit once measured cannot be neasured again. Suppose q G f)y> for the above 
relations to be defined, one needs the indices i, j on which the various command apply to be in V. 
One also needs T to contain the domains of s and t, so that sr and tr are well-defined. This will 
always be the case during the run of a pattern because of condition (D). 

All commands except measurements are deterministic and only modify the quantum part of 
the state. The measurement actions on S are not deterministic, so that these are actually binary 
relations on S, and modify both the quantum and classical parts of the state. The usual convention 
has it that when one does a measurement the resulting state is renormalized and the probabilities 
are associated with the transition. We do not adhere to this convention here, instead we leave the 
states unnormalized. The reason for this choice of convention is that this way, the probability of 
reaching a given state can be read off its norm, and the overall treatment is simpler. As we will 
show later, all the patterns implementing unitary operators will have the same probability for all 
the branches and hence we will not need to carry these probabilities explicitly. 

2 These are actually quadruples of the form (V, W, q, F), unless necessary we will suppress the V and the W. 
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We introduce an additional command called signal shifting: 

q,T X q,T[T(i)+s T /i} 

It consists in shifting the measurement outcome at i by the amount sp- Note that the Z-action leaves 
measurements globally invariant, in the sense that |+ a+7r ), \— a+7T ) = \— a ), \+a}- Thus changing a 
to a + 7r amounts to swapping the outcomes of the measurements, and one has: 

'[Mf] a = Sj°[M?Y (6) 

and signal shifting allows to dispose of the Z action of a measurement, resulting sometimes in 
convenient optimizations of standard forms. 

3.2 Denotational semantics 

Let V be a pattern with computation space V, inputs /, outputs O and command sequence 
A n ...A\. To execute a pattern, one starts with some input state q in fj/, together with the 
empty outcome map 0. The input state q is then tensored with as many |+)s as there are non- 
inputs in V (the N commands), so as to obtain a state in the full space Sjy Then E, M and 
C commands in V are applied in sequence from right to left. We can summarize the situation as 
follows: 

- fio 



<~ HJ® P r£ P r H7® Al...A n r^V\0 

If m is the number of measurements, which is also the number of non outputs, then the run may 
follow 2' m different branches. Each branch is associated with a unique binary string s of length m, 
representing the classical outcomes of the measurements along that branch, and a unique branch 
map A s representing the linear transformation from S)j to Sjo along that branch. This map is 
obtained from the operational semantics via the sequence (g^Tj) with 1 < i < n + 1, such that: 

gi.Ti = q® |+ ...+), 

q n +i = q' ^ o 

A ■ 

and for all i < n : q { , Tj q i+1 , F i+1 . 

Definition 4 A pattern V realizes a map on density matrices p given by p i— ► ^ s At(p)A s . We 
write fPj for the map realized by V. 

Proposition 5 Each pattern realizes a completely positive trace preserving map. 

Proof. Later on we will show that every pattern can be put in a semantically equivalent form 
where all the preparations and entanglements appear first, followed by a sequence of measurements 
and finally local Pauli corrections. Hence branch maps decompose as A s = C S T1 S U, where C s is a 
unitary map over F)o collecting all corrections on outputs, II S is a projection from Sjy to fro rep- 
resenting the particular measurements performed along the branch, and U is a unitary embedding 
from F)i to fjy collecting the branch preparations, and entanglements. Note that U is the same on 
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all branches. Therefore, 

E s 4^s = E s uHclc s u s u 
= J2 s u j nin s u 

= u*(£ a n.)u 
= u^u = i 

where we have used the fact that C s is unitary, Il s is a projection and U is independent of 
the branches and is also unitary. Therefore the map T(p) := ^2 s A s (p)At is a trace-preserving 
completely-positive map (cptp-map), explicitly given as a Kraus decomposition. □ 

Hence the denotational semantics of a pattern is a cptp-map. In our denotational semantics we 
view the pattern as defining a map from the input qubits to the output qubits. We do not explicitly 
represent the result of measuring the final qubits; these may be of interest in some cases. Techniques 
for dealing with classical output explicitly are given by Selinger [Sel04bj and Unruh |Unr05j . 

Definition 6 A pattern is said to be deterministic if it realizes a cptp-map that sends pure states 
to pure states. A pattern is said to be strongly deterministic when branch maps are equal. 

This is equivalent to saying that for a deterministic pattern branch maps are proportional, that 
is to say, for all q G S)i and all Si, S3 € ZJJ - , A Sl (q) and A S2 (q) differ only up to a scalar. For a 
strongly deterministic pattern we have for all Si, S3 G EJ, A Sl = A S2 . 

Proposition 7 If a pattern is strongly deterministic, then it realizes a unitary embedding. 

Proof. Define T to be the map realized by the pattern. We have T = ^ s -Aj-A s . Since the pattern 
in strongly deterministic all the branch maps are the same. Define A to be 2 n / 2 A s , then A must 
be a unitary embedding, because A> A = I. □ 

3.3 Short examples 

For the rest of paper we assume that all the non-input qubits are prepared in the state |+) and 
hence for simplicity we omit the preparation commands Njc 

First we give a quick example of a deterministic pattern that has branches with different proba- 
bilities. Its type is V = {1, 2}, I = O = {1}, and its command sequence is Mr?- Therefore, starting 
with input q, one gets two branches: 

Ma { §(l + e--)g,0[O/2] 
g®|+),0 < 

( i(l-e- fa )g,0[l/2] 

Thus this pattern is indeed deterministic, and implements the identity up to a global phase, and 
yet the two branches have respective probabilities (1 + cosa)/2 and (1 — cosa)/2, which are not 
equal in general and hence this pattern is not strongly deterministic. 

There is an interesting variation on this first example. The pattern of interest, call it T, has 
the same type as above with command sequence X^M^E^. Again, T is deterministic, but not 
strongly deterministic: the branches have different probabilities, as in the preceding example. Now, 
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however, these probabilities may depend on the input. The associated transformation is a cptp-map, 
T(p) := ApA^ + BpB^ with: 

One has A* A + B^B = I, so T is indeed a completely positive and trace-preserving linear map and 
T(\ip){ip\) = (if), ^)|0)(0| and clearly for no unitary U does one have T(p) := UpW. 

For our final example, we return to the pattern 7i, already defined above. Consider the pattern 
with the same qubit space {1,2}, and the same inputs and outputs / = {1}, O = {2}, as 7i, but 
with a shorter command sequence namely M®Ei2- Starting with input q = (a|0) + 6|1))|+), one 
has two computation branches, branching at M®: 

(a|O) + 6|l))|+),0 ^ ^(a|00)+o|01) + 6|10) -6|11)),0 

i((a + 6)|O) + (a-6)|l)),0[O/l] 



M° 



i((a-6)|O) + (a + 6)|l)),0[l/l] 



and since \\a + b\\ 2 + \\a — b\\ 2 = 2(||a|| 2 + ||6|| 2 ), both transitions happen with equal probabilities \. 
Both branches end up with non proportional outputs, so the pattern is not deterministic. However, 
if one applies the local correction X 2 on either of the branches' ends, both outputs will be made to 
coincide. If we choose to let the correction apply to the second branch, we obtain the pattern TC, 
already defined. We have just proved H = U-h, that is to say TL realizes the Hadamard operator. 



3.4 Compositionality of the Denotational Semantics 

With our definitions in place, we will show that the denotational semantics is compositional. 
Theorem 1 For two patterns Vi and V 2 we have {ViPi\ = XP2W1I and Pi®^] = [V2] ® {Vij. 

Proof. Recall that two patterns V\ , Vi may be combined by composition provided V\ has as many 

outputs as V2 has inputs. Suppose this is the case, and suppose further that V\ and V2 respectively 

realize some cptp-maps T\ and T%. We need to show that the composite pattern V2V1 realizes T 2 T\. 
Indeed, the two diagrams representing branches in V\ and V2- 

>')/ •• fio! $)i 2 : - fic;.. 



jO Zl x 11 Si fi Vl x ^ S) 0l x Z^ 1 " 01 Sji 2 x ™ ?)v 2 x Zf Sjo 2 x Z^ 2 "° 2 

can be pasted together, since 0\ = I 2 , and Sjoi — &h- But then, it is enough to notice 1) that 
preparation steps P2 in V2 commute with all actions in V\ since they apply on disjoint sets of qubits, 
and 2) that no action taken in V2 depends on the measurements outcomes in V\. It follows that 
the pasted diagram describes the same branches as does the one associated to the composite V2V1. 

A similar argument applies to the case of a tensor combination, and one has that V2 <8> V\ 
realizes T2 <8) T\. □ 
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If one wanted to give a categorical treatment one can define a category where the objects are 
finite sets representing the input and output qubits and the morphisms are the patterns. This is 
clearly a monoidal category with our tensor operation as the monoidal structure. One can show that 
the denotational semantics gives a monoidal functor into the category of superoperators or into any 
suitably enriched strongly compact closed category |AC04| or dagger category |Sel05aj . It would be 
very interesting to explore exactly what additional categorical structures are required to interpret 
the measurement calculus presented below. Duncan Ross |Dun05| has skectched a polycategorical 
presentation of our measurement calculus. 



4 Universality 

Define the two following patterns on V = {1,2}: 

J (a) := X^M^En (7) 
AZ := E 12 (8) 

with / = {1}, O = {2} in the first pattern, and I = O = {1, 2} in the second. Note that the second 
pattern does have overlapping inputs and outputs. 

Proposition 8 The patterns J (a) and AZ are universal. 

Proof. First, we claim J(a) and AZ respectively realize J(a) and AZ, with: 

We have already seen in our example that J(0) = H implements H = J(0), thus we already know 
this in the particular case where a = 0. The general case follows by the same kind of computation^ 
The case of AZ is obvious. 

Second, we know that these unitaries form a universal set for (gi^C 2 [DKP05j . Therefore, from 
the preceding section, we infer that combining the corresponding patterns will generate patterns 
realizing any unitary in (g> n C 2 . □ 

These patterns are indeed among the simplest possible. As a consequence, in the section devoted 
to examples, we will find that our implementations often have lower space complexity than the 
traditional implementations. 

Remarkably, in our set of generators, one finds a single measurement and a single dependency, 
which occurs in the correction phase of J{a). Clearly one needs at least one measurement, since 
patterns without measurements can only implement unitaries in the Clifford group. It is also true 
that dependencies are needed for universality, but we have to wait for the development of the 
measurement calculus in the next section to give a proof of this fact. 



3 The rest of the paragraph can be omitted without loss of continuity. 

X^Mi a E 12 = X 2 sl Af?P(a)i£i2 = HP(a)i. 



4 Equivalently, this follows from J(a) — HP(a), with P(a) = \ \ ? a ) and: 
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5 The measurement calculus 



We turn to the next important matter of the paper, namely standardization. The idea is quite 
simple. It is enough to provide local pattern-rewrite rules pushing Es to the beginning of the 
pattern and Cs to the end. The crucial point is to justify using the equations as rewrite rules. 



5.1 The equations 

The expressions appearing as commands are all linear operators on Hilbert space. At first glance, 
the appropriate equality between commands is equality as operators. For the deterministic com- 
mands, the equality that we consider is indeed equality as operators. This equality implies equality 
in the denotational semantics. However, for measurement commands one needs a stricter definition 
for equality in order to be able to apply them as rewriting rules. Essentially we have to take into 
the account the effect of different branches that might result from the measurement process. The 
precise definition is below. 

Definition 9 Consider two patterns V and V' we define V = V if and only if for any branch s, 
we have = , where and A£ are the branch map A s defined in Section \3.2[ 

The first set of equations gives the means to propagate local Pauli corrections through the 
entangling operator E^. 

(9) 
(10) 

(11) 
(12) 

These equations are easy to verify and are natural since Eij belongs to the Clifford group, and 
therefore maps under conjugation the Pauli group to itself. Note that, despite the symmetry of the 
Eij operator qua operator, we have to consider all the cases, since the rewrite system defined below 
does not allow one to rewrite E^ to Ej{. If we did allow this the reqrite process could loop forever. 

A second set of equations allows one to push corrections through measurements acting on the 
same qubit. Again there are two cases: 

*[Mf] s A7 = '[Mf] s+r (13) 
\M?\ s Zl = t+r [M«] s (14) 

These equations follow easily from equations (Jl|) and ([5]). They express the fact that the measure- 
ments Mf are closed under conjugation by the Pauli group, very much like equations (|9j) . (fTU|) , (fTT|) 
and (|12p express the fact that the Pauli group is closed under conjugation by the entanglements 
E^. 

Define the following convenient abbreviations: 



EijX- 


V s 7 s TP 
— ^i^j^ij 


EijX- 


V s 7 s TP 
— -X-j^i-^ij 


EijZf 


= ZfEij 


EijZ* 


= Z°Eij 



[M?] s := °[Mf ] s , *[Mf ] := *[Mf ]°, Mf := °[Mf]° 



Mf := MP, Mf := M? 
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Particular cases of the equations above are: 



MfXf = Mf 

MfXf = [MfY = s [Mf] = MfZf 

The first equation, follows from the fact that —0 = 0, so the X action on Mf is trivial; the second 
equation, is because — f is equal f + vr modulo 2ir, and therefore the X and Z actions coincide on 
Mf . So we obtain the following: 

l [Mf] s = *[Mf] (15) 
*[MJT = ^[M?] (16) 

which we will use later to prove that patterns with measurements of the form M x and M y may 
only realize unitaries in the Clifford group. 

5.2 The rewrite rules 

We now define a set of rewrite rules, obtained by orienting the equations abov^f): 

EijXf =>• XfZjEij EX 

EijXj => XjZfEij EX 

EijZf =>• ZfEij EZ 

EijZ? =>■ Z?Eij EZ 

*[Mf] s A7 =► *[Mf] s+r MX 

*[Mf] s Z[ =► r+ *[Mf] s MZ 

to which we need to add the free commutation rules, obtained when commands operate on disjoint 
sets of qubits: 

EijA^ => A-X^ij where A is not an entanglement 
ArXf =>■ Xf At where A is not a correction 
ArZf =>■ ZfAf where ^4 is not a correction 

k i « fe 

where A; represent the qubits acted upon by command A, and are supposed to be distinct from i 
and j. Clearly these rules could be reversed since they hold as equations but we are orienting them 
this way in order to obtain termination. 

Condition (D) is easily seen to be preserved under rewriting. 

Under rewriting, the computation space, inputs and outputs remain the same, and so do the 
entanglement commands. Measurements might be modified, but there is still the same number 
of them, and they still act on the same qubits. The only induced modifications concern local 
corrections and dependencies. If there was no dependency at the start, none will be created in the 
rewriting process. 

In order to obtain rewrite rules, it was essential that the entangling command (AZ) belongs 
to the normalizer of the Pauli group. The point is that the Pauli operators are the correction 
operators and they can be dependent, thus we can commute the entangling commands to the 
beginning without inheriting any dependency. Therefore the entanglement resource can indeed be 
prepared at the outset of the computation. 

5 Recall that patterns are executed from right to left. 
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5.3 Standardization 



Write V V , respectively V V, if both patterns have the same type, and one obtains the 
command sequence of V' from the command sequence of V by applying one, respectively any 
number, of the rewrite rules of the previous section. We say that V is standard if for no V' , V =>■ V' 
and the procedure of writing a pattern to standard form is called standardization^. 

One of the most important results about the rewrite system is that it has the desirable properties 
of determinacy (confluence) and termination (standardization). In other words, we will show that 
for all V , there exists a unique standard V' , such that V V 1 . It is, of course, crucial that the 
standardization process leaves the semantics of patterns invariant. This is the subject of the next 
simple, but important, proposition, 

Proposition 10 Whenever V V , \V\ = {V'j. 

Proof. It is enough to prove it when V =>■ V . The hrst group of rewrites has been proved to be 
sound in the preceding subsections, while the free commutation rules are obviously sound. □ 

We now begin the main proof of this section. First, we prove termination. 

Theorem 2 (Termination) All rewriting sequences beginning with a pattern V terminate after 
finitely many steps. For our rewrite system, this implies that for all V there exist finitely many V 1 
such that V V where the V' are standard. 

Proof. Suppose V has command sequence A n . . . A\\ so the number of commands is n. Let 
e < n be the number of E commands in V . As we have noted earlier, this number is invariant 
under =x Moreover E commands in V can be ordered by increasing depth, read from right to left, 
and this order, written <e, is also invariant, since EE commutations are forbidden explicitly in 
the free commutation rules. 

Define the following depth function d on E and C commands in V: 



Define further the following sequence of length e, dEiV){i) is the depth of the E'-command of rank 
i according to <e- By construction this sequence is strictly increasing. Finally, we define the 
measure m(V) := {dE{V),dc{V)) with: 



We claim the measure we just defined decreases lexicographically under rewriting, in other words 
V => V implies m(V) > m(V'), where < is the lexicographic ordering on N e+1 . 

To clarify these definitions, consider the following example. Suppose V's command sequence is 
of the form EXZE, then e = 2, cIe(V) = (1,4), and m(V) = (1,4,3). For the command sequence 
EEX we get that e = 2, d,E(V) = (2,3) and m(V) = (2,3,2). Now, if one considers the rewrite 
EEX =>■ EXZE, the measure of the left hand side is (2, 3, 2), while the measure of the right hand 
side, as said, is (1,4,3), and indeed (2,3,2) > (1,4,3). Intuitively the reason is clear: the Cs are 
being pushed to the left, thus decreasing the depths of Es, and concomitantly, the value of dg- 

6 We use the word "standardization" instead of the more usual "normalization" in order not to cause terminological 
confusion with the physicists' notion of normalization. 




dc(V) 
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Let us now consider all cases starting with an EC rewrite. Suppose the E command under 
rewrite has depth d and rank i in the order <e- Then all Es of smaller rank have same depth in 
the right hand side, while E has now depth d — 1 and still rank i. So the right hand side has a 
strictly smaller measure. Note that when C = X, because of the creation of a Z (see the example 
above), the last element of m('P) may increase, and for the same reason all elements of index j > i 
in dEiV) may increase. This is why we are working with a lexicographical ordering. 

Suppose now one does an MC rewrite, then dc(V) strictly decreases, since one correction is 
absorbed, while all E commands have equal or smaller depths. Again the measure strictly decreases. 

Next, suppose one does an EA rewrite, and the E command under rewrite has depth d and rank 
i. Then it has depth d— 1 in the right hand side, and all other E commands have invariant depths, 
since we forbade the case when A is itself an E. It follows that the measure strictly decreases. 

Finally, upon an AC rewrite, all E commands have invariant depth, except possibly one which 
has smaller depth in the case A = E, and dc(V) decreases strictly because we forbade the case 
where A = C. Again the claim follows. 

So all rewrites decrease our ordinal measure, and therefore all sequences of rewrites are finite, 
and since the system is finitely branching (there are no more than n possible single step rewrites 
on a given sequence of length n), we get the statement of the theorem. 

The final statement of the theorem follows from the fact that we have finitely many rules so the 
system is finitely branching. In any finitely branching rewrite system with the property that every 
rewrite sequence terminates, it is clearly true that there can be only finitely many standard forms. 
□ 

The next theorem establishes the important determinacy property and furthermore shows that 
the standard patterns have a certain canonical form which we call the NEMC form. The precise 
definition is: 

Definition 11 A pattern has a NEMC form if its commands occur in the order of Ns first, then 
Es , then Ms, and finally Cs. 

We will usually just say "EMC" form since we can assume that all the auxiliary qubits are prepared 
in the |+) state we usually just elide these N commands. 

Theorem 3 (Confluence) For all V , there exists a unique standard V ' , such that V =>* V , and 
V' is in EMC form. 

Proof. Since the rewriting system is terminating, confluence follows from local confluence by 
Newman's lemma, see, for example, [Bar 84]. The uniqueness of the standard is form an immediate 
consequence. 

We look for critical pairs, that is occurrences of three successive commands where two rules can 
be applied simultaneously. One finds that there are only five types of critical pairs, of these the 
three involve the N command, these are of the form: NMC, NEC and NEM; and the remaining 
two are: EijMkCk with i, j and k all distinct, EijM^Ci with k and I distinct. In all cases local 
confluence is easily verified. 

Suppose now V 1 does not satisfy the EMC form conditions. Then, either there is a pattern EA 
with A not of type E, or there is a pattern AC with A not of type C. In the former case, E and 

7 This means that whenever two rewrite rules can be applied to a term t yielding t\ and fe, one can rewrite both 
ti and t2 to a common third term £3, possibly in many steps. 
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A must operate on overlapping qubits, else one may apply a free commutation rule, and A may 
not be a C since in this case one may apply an EC rewrite. The only remaining case is when A 
is of type M, overlapping E's qubits, but this is what condition (Dl) forbids, and since (Dl) is 
preserved under rewriting, this contradicts the assumption. The latter case is even simpler. □ 

We have shown that under rewriting any pattern can be put in EMC form, which is what we 
wanted. We actually proved more, namely that the standard form obtained is unique. However, 
one has to be a bit careful about the significance of this additional piece of information. Note first 
that uniqueness is obtained because we dropped the CC and EE free commutations, thus having a 
rigid notion of command sequence. One cannot put them back as rewrite rules, since they obviously 
ruin termination and uniqueness of standard forms. 

A reasonable thing to do, would be to take this set of equations as generating an equivalence 
relation on command sequences, call it =, and hope to strengthen the results obtained so far, by 
proving that all reachable standard forms are equivalent. 

But this is too naive a strategy, since E\<iX 1 X 2 = E\2X 2 X 1 , and: 

I? ys v~t * vs rys vt ryt TP 

VS yt 17S Vt TP 

obtaining an expression which is not symmetric in 1 and 2. To conclude, one has to extend = 
to include the additional equivalence XfZ\ = Z^Xf, which fortunately is sound since these two 
operators are equal up to a global phase. Thus, these are all equivalent in our semantics of patterns. 
We summarize this discussion as follows. 

Definition 12 We define an equivalence relation = on patterns by taking all the rewrite rules as 
equations and adding the equation X\Z\ = Z\X\ and generating the smallest equivalence relation. 

With this definition we can state the following proposition. 

Proposition 13 All patterns that are equivalent by = are equal in the denotational semantics. 

This = relation preserves both the type (the (V, /, O) triple) and the underlying entanglement graph. 
So clearly semantic equality does not entail equality up to =. In fact, by composing teleportation 
patterns one obtains infinitely many patterns for the identity which are all different up to =. One 
may wonder whether two patterns with same semantics, type and underlying entanglement graph 
are necessarily equal up to =. This is not true either. One has J(a)J(0)J((3) = J (a + (3) = 
J(P)J(0)J(a) (where J(a) is defined in Section H]), and this readily gives a counter-example. 
We can now formally describe a simple standardization algorithm. 

Algorithm 1 Input: A pattern V on |V| = N qubits with command sequence Am ••• A\. 
Output: An equivalent pattern V 1 in NEMC form. 

1. Commute all the preparation commands (new qubits) to the right side. 

2. Commute all the correction commands to the left side using the EC and MC rewriting rules. 

3. Commute all the entanglement commands to the right side after the preparation commands. 

Note that since each qubit can be entangled with at most N — 1 other qubits, and can be 
measured or corrected only once, we have 0(N 2 ) entanglement commands and O(N) measurement 
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commands. According to the definiteness condition, no command acts on a qubit not yet prepared, 
hence the first step of the above algorithm is based on trivial commuting rules; the same is true 
for the last step as no entanglement command can act on a qubit that has been measured. Both 
steps can be done in 0{N 2 ). The real complexity of the algorithm comes from the second step 
and the EX commuting rule. In the worst case scenario, commuting an X correction to the left 
might create 0(N 2 ) other Z corrections, each of which has to be commuted to the left themselves. 
Thus one can have at most 0(N 3 ) new corrections, each of which has to be commuted past 0(N 2 ) 
measurement or entanglement commands. Therefore the second step, and hence the algorithm, has 
a worst case complexity of 0(N 5 ). 

We conclude this subsection by emphasizing the importance of the EMC form. Since the 
entanglement can always be done first, we can always derive the entanglement resource needed for 
the whole computation right at the beginning. After that only local operations will be performed. 
This will separate the analysis of entanglement resource requirements from the classical control. 
Furthermore, this makes it possible to extract the maximal parallelism for the execution of the 
pattern since the necessary dependecies are explicitly expressed, see the example in section [6] for 
further discussion. Finally, the EMC form provides us with tools to prove general theorems about 
patterns, such as the fact that they always compute cptp-maps and the expressiveness theorems of 
section [7J 



5.4 Signal shifting 

One can extend the calculus to include the signal shifting command Sj. This allows one to dispose 
of dependencies induced by the Z-action, and obtain sometimes standard patterns with smaller 
computational depth complexity, as we will see in the next section which is devoted to examples. 



[M?Y Sj[M? 

ysl > sjxf 



where s[t/si] denotes the substitution of Sj with t in s, s, t being signals. Note that when we write 
a t explicitly on the upper left of an M, we mean that t ^ 0. The first additional rewrite rule was 
already introduced as equation ([6]), while the other ones merely propagate the signal shift. Clearly 
one can dispose of Sj when it hits the end of the pattern command sequence. We will refer to this 
new set of rules as =>s- Note that we always apply first the standardization rules and then signal 
shifting, hence we do not need any commutation rule for E and S commands. 

It is important to note that both theorem [2] and [3] still hold for this extended rewriting system. 
In order to prove termination one can start with the EMC form and then adapt the proof of 
Theorem [2] by defining a depth function for a signal shift similar to the depth of a correction 
command. As with the correction, signal shifts can also be commuted to the left hand side of a 
command sequence. Now our measure can be modified to account for the new signal shifting terms 
and shown to be decreasing under each step of signal shifting. Confluence can be also proved from 
local confluence using again Newman's Lemma |Bar84| . One typical critical pair is t [Mj']Sf where 
i appears in the domain of signal t and hence the signal shifting command Sf will have an effect on 
the measurement. Now there are two possible ways to rewrite this pair, first, commute the signal 
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shifting command and then replace the left signal of the measurement with its own signal shifting 
command: 

*[M?] S? => St t+s [Mf) 
=> S?S° +t M? 

The other way is to first replace the left signal of the measurement and then commute the signal 
shifting command: 

*[Af?]5? => SjMfSf 
S j S? Mf 

Now one more step of rewriting on the last equation will give us the same result for both choices. 

S^SfMf => SfS^Mf 
All other critical terms can be dealt with similarly. 



6 Examples 

In this section we develop some examples illustrating pattern composition, pattern standardization, 
and signal shifting. We compare our implementations with the implementations given in the refer- 
ence paper [RBB03] . To combine patterns one needs to rename their qubits as we already noted. 
We use the following concrete notation: if V is a pattern over {1, . . . , n}, and / is an injection, 
we write V(f(l),...,f(n)) for the same pattern with qubits renamed according to /. We also 
write V 2 o V\ for pattern composition, in order to make it more readable. Finally we define the 
computational depth complexity to be the number of measurement rounds plus one final correction 
round. More details on depth complexity, specially on the preparation depth, i.e. depth of the 
entanglement commands, can be found in |BK06] . 



Teleportation. 



Consider the composite pattern J((3){2, 3)o l 7(a)(l, 2) with computation space {1, 2, 3}, inputs {1}, 
and outputs {3}. We run our standardization procedure so as to obtain an equivalent standard 
pattern: 



rS 2 

L 3 



X? 





E23X2 1 


M~ P X^ 





M^ a E 12 



^23-^12 



Let us call the pattern just obtained J(ot, 13). If we take as a special case a = (3 = 0, we get: 

X S 3 2 Z S 3 1 M$M?E 23 E 12 

and since we know that J7(0) implements H and H 2 = I, we conclude that this pattern implements 
the identity, or in other words it teleports qubit 1 to qubit 3. As it happens, this pattern obtained 
by self-composition, is the same as the one given in the reference paper |RBB03t p. 14]. 
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x-rotation. 



Here is the reference implementation of an x-rotation [RBB031 p. 17], R x { 



XI 2 Z* 1 [M^ a ]^MfE 23 E 12 

with type {1,2,3}, {1}, and {3}. There is a natural question which one might call the recognition 
problem, namely how does one know this is implementing R x (a) ? Of course there is the brute 
force answer to that, which we applied to compute our simpler patterns, and which consists in 
computing down all the four possible branches generated by the measurements at qubits 1 and 2. 
Another possibility is to use the stabilizer formalism as explained in the reference paper [RBB03] . 
Yet another possibility is to use pattern composition, as we did before, and this is what we are 
going to do. 

We know that R x {a) = J{a)H up to a global phase, hence the composite pattern l 7(a)(2, 3) o 
"H(l,2) implements R x {a). Now we may standardize it: 



J(a)(2,3)oW(l,2) 



'EX 



'MX 



A 3 Z 3 



E 23 X 2 1 



AC 



MfE 12 
M£E 23 Ei 2 



X^Z* 1 [M^ a ] Sl MfE 23 E 
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obtaining exactly the implementation above. Since our calculus preserves the semantics, we deduce 
that the implementation is correct. 

z-rotation. 

Now, we have a method here for synthesizing further implementations. Let us replay it with another 
rotation R z {a). Again we know that R z (a) = HR x {a)H , and we already know how to implement 
both components H and R x (a). 

So we start with the pattern W(4, 5) o TZ x (a)(2, 3, 4) o H(l, 2) and standardize it: 



ft(4,5)o^.(a)(2,3,4)oH(l,2) 



W(4, 5)X| 3 Zf [M 3 a ] 1+S2 Mf £ 3 4 



E23X0 1 



H{4, 5)Xf Zf [M 3 a ] 1+S2 Mf 
W(4, 5)X% 3 Zf [M 3 a ] 1+S2 Z^ 



X| 4 M£ 



ft(4,5)Xf Zf [M 3 a ] 1+ * 2 Zf MZMfEuM 



E45X4 3 



M£Xf 



M?E 12 
M?E 123 



Mf £"1234 



Zf Sl [M 3 Q ] 1+S2 Mf Mf E1234 

MfAff£?i2345 



Zf Sl [M; 



aT|l+S2 



^[Mi) 1+S2 M^MfE 12M5 



'EX 

'EZ 

'MX 

'MZ 

'EX 

'MX 

'MZ 



X S4 Z s 3 S2 [ M |] S3 « [ M a] 1+S2 Mf Mf £12345 



To aid reading E 23 E\ 2 is shortened to -E123, E\ 2 E 23 E 3 ^ to £1234, and *[M"] 1+S is used as shorthand 
for *[M" Q ] S . 

Here for the first time, we see MZ rewritings, inducing the Z-action on measurements. The 
resulting standardized pattern can therefore be rewritten further using the extended calculus: 

X^Z 3 ^^ [Mf] 8 ^ [Mgaf+^Mf Mf £ 123 45 
X s 2 +S4 z' x+a *M% [M 3 Q ] 1+52 Mf Mf £12345 



21 



obtaining the pattern given in the reference paper [RBB031 p. 5]. 

However, just as in the case of the R x rotation, we also have R z (a) = HJ(a) up to a global 
phase, hence the pattern TC(2, 3)J'(a)(l, 2) also implements R z (a), and we may standardize it: 



W(2,3)o J(a)(l,2) 



X| 2 Mf 

A 3 Z 3 



^mx X s 3 2 Z^M%M^ a E 123 



E23X0 1 



M{ a E 12 
M{ a E 123 



obtaining a 3 qubit standard pattern for the z-rotation, which is simpler than the preceding one, 
because it is based on the J (a) generators. Since the z-rotation R z (a) is the same as the phase 
operator: 

= (J x 

up to a global phase, we also obtain with the same pattern an implementation of the phase oper- 
ator. In particular, if a = f, using the extended calculus, we get the following pattern for -P(f): 
X^ 2 Z^ 1+1 M^MfE 123 . 



General rotation. 

The realization of a general rotation based on the Euler decomposition of rotations as R x (^)R z {(3)R x (a) 
would results in a 7 qubit pattern. We get a 5 qubit implementation based on the J(a) decompo- 
sition [DKP05| : 

7) = J(0)J(-a)J(-P)J(-i) 

(The parameter angles are inverted to make the computation below more readable.) The extended 
standardization procedure yields: 

J(0)(A, 5)J(-a)(3, 4) J(-/3)(2, 3)J(- j)(l, 2) = 
Xl 4 M2E 45 X^M^E 34 Xl 2 Mj\ E 23 X^ \ M?E 12 

X^MlE^X^MiE^xf 

x s 5 4 m2e A5 xi 3 m% 



x s 5 *m2e 45 x s / 



x s 5 4 m2 



M^X* 1 



-C/34A3 Z 3 



M^X^Z* 1 



E/^XTZl 2 



X^Z^ 3S2 



MlX^Zl* 



Z s ^MlE 123 
\M^MlE X23 



f] 

Z? [M^M^E 1234 
W§] S2 [M^M^E 1234 
Z^[M«] S2 [M^M^E 12U5 
M^[M^M^E 12345 



> EX 
> MX 
> EXZ 
>MXZ 
> EXZ 
> MXZ 
'S 



xi 2+S4 z s 5 1+S3 m2 [m 3 °] S2 [m^] si m^e 12345 



CNOT (AX). 

This is our first example with two inputs and two outputs. We use here the trivial pattern X with 
computation space {1}, inputs {1}, outputs {1}, and empty command sequence, which implements 
the identity over f)\. 
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One has AX = (I <g) H) f\Z(I (g) i7), so we get a pattern using 4 qubits over {1, 2, 3, 4}, with 
inputs {1,2}, and outputs {1,4}, where one notices that inputs and outputs intersect on the control 
qubit {1}: 

(1(1) ® ((3, 4)) AZ(l, 3)(T(1) ® ((2, 3)) = X s 4 3 Mf E U E 13 X S 3 2 M$ E 23 

By standardizing: 

£ 13 Xg 2 | Mf E 23 

-E'34^3 2 Mjf E13E23 =^EX 
Mf Xg 2 | Mf ^13^23^34 ^AfX 

XfZfZfMfMf E 13 E 23 E U 

Note that, in this case, we are not using the -E1234 abbreviation, because the underlying struc- 
ture of entanglement is not a chain. This pattern was already described in Aliferis and Leung's 
paper |AL04| . In their original presentation the authors actually use an explicit identity pattern (us- 
ing the teleportation pattern J(0, 0) presented above), but we know from the careful presentation 
of composition that this is not necessary. 



X S 4 3 M*E U 
X S A Z Z{ 2 Mf 

VS3 ?S2 7S2 
^4 -^4 ^1 



GHZ. 

We present now a family of patterns preparing the GHZ entangled states |0 . . . 0) + |1 . . . 1). One 
has: 

GHZ(n) = (H n AZ n „ ln ...H 2 AZ 12 ) | + ...+) 

and by combining the patterns for AZ and H, we obtain a pattern with computation space 
{1, 2, 2', ... , n, n'}, no inputs, outputs {1, 2', . . . , n'}, and the following command sequence: 

X^?M*E nn tE( n _iy n . . . XyM%E 22 iE\ 2 

With this form, the only way to run the pattern is to execute all commands in sequence. The 
situation changes completely, when we bring the pattern to extended standard form: 



X s jMZE nn ,E {n _ iyn . . . X s 3 fM$E 33 , 


E 2 < 3 Xy 


M 2 E 22 iE\ 2 


X n' ''X 2 ? : M£E nn > -E( n _i)/ n • • • X°? 




M 2 E 33 i E 2 / 3 E 22 i E\ 2 



X s jX s 2 ?MZE nn ,E in _ lYn . . . X 3 r 2 [Mf ]Mf E 33 ,E 2 , 3 E 22 ,E 12 ^ 
X s n 7... X s 3 fX s 2 f s ^ [M*\ ... 52 [M$]MZE nn ,E {n _ iyn . . . E 33 ,E 2 , 3 E 22 ,E 12 ^ s 
x s, +S3+ ... +Sn X 82+s ax g M x . . . MlMlE nn ,E {n _ xyn . . . E 33 ,E 2 , 3 E 22 ,E l2 

All measurements are now independent of each other, it is therefore possible after the entanglement 
phase, to do all of them in one round, and in a subsequent round to do all local corrections. In 
other words, the obtained pattern has constant computational depth complexity 2. 



Controlled-C/. 

This final example presents another instance where standardization obtains a low computational 
depth complexity, the proof of this fact can be found in [BK06J . For any 1-qubit unitary U, one 
has the following decomposition of All in terms of the generators J (a) [DKP 05J : 

7 7T 7T 7 — 7T — 6 — /3 —0 + 6 — TT 

AC/12 = Ji Jf ' J 2 J 2 J 2 2 J 2 2 J 2 AZ12 J 2 2 J 2 J 2 2 J$AZ 12 J 2 2 
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with a' 



a + 



/3+7+<S 



By translating each J operator to its corresponding pattern, we get the 



following wild pattern for f\U : 



xy>M%EBcX% A MZ*E AB X' k i M9E jk X a j t M7f l -*E ij 
X? h Ml E hl X s h s M^ E gh X s g f M° f E fg E A f X s f M~ 5 E ef 

7T + 5 + j3 — 5 + TT 

X s *M*E de X s /M c 2 E cd X s *MlE hc E Ah X s h *M a 2 E ab 



In order to run the wild form of the pattern one needs to follow the pattern commands in sequence. 
It is easy to verify that, because of the dependent corrections, one needs at least 12 rounds to 
complete the execution of the pattern. The situation changes completely after extended standard- 



ization: 



7 Si + Sg + Se+Sc+Sa v $j +Sh + S f +Sd+Sb V SB rySA+Se+S 



7 tt+S+13 



s »M?M a 2 



J k yl k 

M^M^ a ' Mj [M- 

M° f [M~^} Sd+Sb [M~ 
Ebc Eab Ej & Eij Eu E g hEf g EA/E e f E de E c dEb c E a b EAb 

Now the order between measurements is relaxed, as one sees in Figure [21 which describes the depen- 
dency structure of the standard pattern above. Specifically, all measurements can be completed in 
7 rounds. This is just one example of how standardization lowers computational depth complexity, 
and reveals inherent parallelism in a pattern. 




Figure 2: The dependency graph for the standard AU pattern. 



7 The no dependency theorems 

From standardization we can also infer results related to dependencies. We start with a simple 
observation which is a direct consequence of standardization. 

Lemma 14 Let V be a pattern implementing some cptp-maps T, and suppose V 's command se- 
quence has measurements only of the M x and M y kind, then U has a standard implementation, 
having only independent measurements, all being of the M x and M v kind (therefore of computa- 
tional depth complexity at most 2). 
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Proof. Write V' for the standard pattern associated to V . By equations (fT5j) and (fTBj) . the X- 
actions can be eliminated from V' , and then Z-actions can be eliminated by using the extended 
calculus. The final pattern still implements T, has no longer any dependent measurements, and 
has therefore computational depth complexity at most 2. □ 

Theorem 4 Let U be a unitary operator, then U is in the Clifford group iff there exists a pattern 
V implementing U , having measurements only of the M x and M y kind. 

Proof. The "only if" direction is easy, since we have seen in the example section, standard patterns 
for AX, H and which had only independent M x and M y measurements. Hence any Clifford 

operator can be implemented by a combination of these patterns. By the lemma above, we know 
we can actually choose these patterns to be standard. 

For the "if" direction, we prove that U belongs to the normalizer of the Pauli group, and hence 
by definition to the Clifford group. In order to do so we use the standard form of V written as 
V' = C-piM-piE-pi which still implements U, and has only M x and M y measurements. Recall that, 
because of equations (fT5|) and (fT6l) . these measurements are independent. 

Let i be an input qubit, and consider the pattern V" = V'Ci, where Cj is either Xj or Z{. 
Clearly V" implements UCi. First, one has: 

C-pi M'pi E'piCi C'ptM'prC 1 ' Epi 

for some non-dependent sequence of corrections C", which, up to free commutations can be written 
uniquely as C'qC" , where C'q applies on output qubits, and therefore commutes to M-pi, and C" 
applies on non-output qubits (which are therefore all measured in Mp/). So, by commuting C' Q 
both through Mp/ and C-p> (up to a global phase), one gets: 

Using equations (fT5j) . (fTBT) . and the extended calculus to eliminate the remaining Z-actions, one 
gets: 

M<p,C" ^mc,s SM v 

for some product S = Y[{j^j} &j °f constant shifts applying to some subset J of the non-output 
qubits. So: 

CqC-p' M-p/C Epi s C qCj>i S M-pi Epi 

CqCqCp' M-pi E-pi 

where Cq is a further constant correction obtained by signal shifting C-pi with 5. This proves that 
V" also implements C'qCqU, and therefore UCi = C'qCqU which completes the proof, since C'qC'q 
is a non dependent correction. □ 

The "only if" part of this theorem already appears in previous work [RBB03, p. 18]. The "if" 
part can be construed as an internalization of the argument implicit in the proof of Gottesman-Knill 
theorem [NC001 p.464]. 

We can further prove that dependencies are crucial for the universality of the model. Observe 
first that if a pattern has no measurements, and hence no dependencies, then it follows from (D2) 
that V = O, i.e., all qubits are outputs. Therefore computation steps involve only X, Z and 

8 Here we have used the trivial equations Z°- +1 = ZiZI and = XiXf 
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AZ, and it is not surprising that they compute a unitary which is in the Clifford group. The 
general argument essentially consists in showing that when there are measurements, but still no 
dependencies, then the measurements are playing no part in the result. 

Theorem 5 Let V be a pattern implementing some unitary U, and suppose V 's command sequence 
doesn't have any dependencies, then U is in the Clifford group. 

Proof. Write V' for the standard pattern associated to V . Since rewriting is sound, V still 
implements U, and since rewriting never creates any dependency, it still has no dependencies. In 
particular, the corrections one finds at the end of V' , call them C, bear no dependencies. Erasing 
them off V 1 , results in a pattern V" which is still standard, still deterministic, and implementing 
U' := C^U. 

Now how does the pattern V" run on some input <f> ? First 4"S> |+ . . . +) goes by the entanglement 
phase to some ip G Sjy, and is then subjected to a sequence of independent 1-qubit measurements. 
Pick a basis B spanning the Hilbert space generated by the non-output qubits Sjy^O an d associated 
to this sequence of measurements. 

Since £jy = fto <8> £)V\0 and fiv^o = ©(fees [</>&]) where [4>b] is the linear subspace generated by 
<f>b, by distributivity, ip uniquely decomposes as: 

^ = E«feeS X b ® <fo 

where ranges over B, and x& G f)o- Now since V" is deterministic, there exists an x, and scalars 
Xb such that Xb = \x. Therefore ip can be written x <8> ip' , for some ip' . It follows in particular 
that the output of the computation will still be x (up to a scalar), no matter what the actual 
measurements are. One can therefore choose them to be all of the M x kind, and by the preceding 
theorem U' is in the Clifford group, and so is U = CU', since C is a Pauli operator. □ 

From this section, we conclude in particular that any universal set of patterns has to include 
dependencies (by the preceding theorem), and also needs to use measurements M a where a / 
modulo 7£ (by the theorem before). This is indeed the case for the universal set J{a) and AZ. 

8 Other Models 

There are several other approaches to measurement-based computation as we have mentioned in 
the introduction. However, it is only for the one-way model that the importance of having all 
the entanglement in front has been emphasized. For example, Gottesman and Chuang describe 
computing with teleportation in the setting of the circuit model and hence the computation is very 
sequential |GC99j . What we will do is to give a general treatment of a variety of measurement- 
based models - including some that appear here for the first time - in the setting of our calculus. 
More precisely we would like to know other potential definitions for commands N, E, M and C 
that lead to a model that still satisfies the properties of: (i) being closed under composition; (ii) 
universality and (iii) standardization. 

Moreover we are interested in obtaining a compositional embedding of these models into a single 
one-qubit measurement-based model. The teleportation model can indeed be embedded into the 
one-way model. There is, however, a new model, the Pauli model - formally defined here for the 
first time - which is motivated by considerations of fault tolerance [RAB04, DK05b. DKOSQ6]. The 
Pauli model can be embedded into a slight generalization of the one-way model called the phase 
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model; also given here for the first time. The one-way model will trivially embed in the phase 
model so by composition all the measurement-based models will embed in the phase model. We 
could have done everything ab initio in terms of the phase model but this would have made much 
of the presentation unnecessarily complicated at the outset. 

We recall the remark from the introduction that these embeddings have three advantages: first, 
we get a workable syntax for handling the dependencies of operators on previous measurement 
outcomes, second, one can use these embeddings to transfer the measurement calculus previously 
developed for the one-way model to obtain a calculus for the new model including, of course, a 
standardization procedure that we get automatically; lastly, one can embed the patterns from the 
phase model into the new models and vice versa. In essence, these compositional embeddings will 
allow us to exhibit the phase model as being a core calculus for measurement-based computation. 
However different models are interesting from the point of view of implementation issues like fault- 
tolerance and ease of preparation of entanglement resources. Our embeddings allow one to move 
easily between these models and to concentrate on the one-way model for designing algorithms and 
proving general theorems. 

This section has been structured into several subsections, one for each model and its embedding. 



8.1 Phase Model 

In the one-way model the auxiliary qubits are initialized to be in the |+) state. We extend the 
one-way model to allow the auxiliary qubits to be in a more general state. We define the extended 
preparation command iV" to be the preparation of the auxiliary qubit i in the state |+«). We 
also add a new correction command Zf, called a phase correction to guarantee that we can obtain 
determinate patterns. The dependent phase correction is written as Z^ ,s with Z®' = I and 

Z"' 1 = ^ J* Q ^ • Under conjugation, the phase correction, defines a new action over measurement: 

(Z /3 ' s ^M a zP' s = M a ~ s/3 

and since the measurement is destructive, it simplifies to M"Z? ,S = M^~ s ^ . This action does 
not commute with Pauli actions and hence one cannot write a compact notation for dependent 
measurement, as we did before, and the computation of angle dependencies is a bit more compli- 
cated. Thereafter, a measurement preceded by a sequence of corrections on the same qubit will be 
called a dependent measurement. Note that, by the absorption equations, this indeed can be seen 
as a measurement, where the angle depends on the outcomes of some other measurements made 
beforehand. 

To complete the extended calculus it remains to define the new rewrite rules: 



EijZ?' 8 -- 






EP 


MfXf -- 


* M t 


-l) s a 


MX 


M?Z? -- 


* K 


— S7V 


MZ 


MfZ^ s -- 


* K 


-s/3 


MP 



The above rules together with the rewriting rules of the one-way model described in Section (5J lead 
to a standardization procedure for the model. It is trivial that the one-way model is a fragment of 
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this generalized model and hence universality immediately follows. It is also easy to check that the 
model is closed under composition and all the semantical properties of the one-way model can be 
extended to this general model as well. 

The choice of extended preparations and its concomitant phase correction is actually quite 
delicate. One wishes to keep the standardizability of the calculus which constrains what can be 
added but one also wishes to have determinate patterns which forces us to put in appropriate 
corrections. The phase model is only a slight extension of the original one-way model, but it allows 
a discussion of the next model which is of great physical interest. 



8.1.1 Pauli model 

An interesting fragment of the phase model is defined by restricting the angles of measurements to 
{0, |, 7r, — ^} i.e. Pauli measurements and the angles of preparation to and \. Also the correction 

7T 

commands are restricted to Pauli corrections X, Z and Phase correction Z». One readily sees that 
the subset of angles is closed under the actions of the corrections and hence the Pauli model is 
closed under composition. 

Proposition 15 The Pauli model is approximately universal. 



Proof. We know that the set consisting of J(0) (which is H), J(j), and AZ is approximately 
universal. Hence, to prove the approximate universality of Pauli model, it is enough to exhibit a 
pattern in the Pauli model for each of these three unitaries. We saw before that J(0) and AZ are 
computed by the following 2-qubit patterns: 



J(0) 
AZ 



X^M^E 12 
E\2 



where both belong also to the Pauli model. The pattern for J| in the one-way model is expressed 
as follows: 

J{\) := X^M^En^ 
= X s 2 l M*E 12 z} 

The above forms do not fit in the Pauli model, since the first one uses a measurement with an angle 
j and the second uses Z*. However by teleporting the input qubit and then applying the Z^ and 
finally running the standardization procedure we obtain the following pattern in the Pauli model 
for J(f ): 

X^M^EuZf 
= X S /M°E 3A _ 



ZfXpMgMfEnEha 



XI 3 M$E U Z° 2 Z 3 2 ' 2 X^ M$M°E 12 E 23 



Xl 3 M^E u Z^ 2 Z^ 



X S ^M^M^E 12 E 23 



X| 3+S2 Zf M 3 



(-i) s i S2 f 



M$M?E 12 E 23 E U N* 



Approximate universality for the Pauli model is now immediate. □ 

Note that we cannot really expect universality (as we had for the phase model) because the 
angles are restricted to a discrete set. On the other hand it is precisely this restriction that makes 
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the Pauli model interesting from the point of view of implementation. The other particular interest 
behind this model, apart from its simple structure, is based on the existence of a novel fault tolerant 
technique for computing within this framework |BK05l IRAB041 IDKOS06] . 

8.2 Teleportation 

Another class of measurement-based models - older, in fact, than the one-way model - uses 2- 
qubit measurements. These are collectively referred to as teleportation models |Leu04] . Several 
papers that are concerned with the relation and possible unification of these models |CLN05t 
IAL044 IJP05] have already appeared. One aspect of these models that stands in the way of a 
complete understanding of this relation, is that, whereas in the one-way model one has a clearly 
identified class of measurements, there is less agreement concerning which measurements are allowed 
in teleportation models. 

We propose here to take as our class of 2-qubit measurements a family obtained as the conjugate 
under the operator AZ of tensors of 1-qubit measurements. We show that the resulting teleportation 
model is universal. Moreover, almost by construction, it embeds into the one-way model, and thus 
exposes completely the relation between the two models. 

Before embarking on the specifics of our family of 2-qubit measurements, we remark that the 
situation commented above is more general: 

Lemma 16 Let A be an orthonormal basis in <S> n C 2 , with associated n-qubit measurement M A , 
and Ai with i = 1, . . . ,n be orthonormal bases in C 2 , with associated 1-qubit measurements M^ 1 . 
Then there exists a unique (up to a permutation) n-qubit unitary operator U such that: 

Mtn = U^n^M^Ui.^ 

Proof. Take U to map ®iAi to A. □ 

This simple lemma says that general n-qubit measurements can always be seen as conjugated 
1-qubit measurements, provided one uses the appropriate unitary to do so. As an example consider 
the orthogonal graph basis Q = AZ\ 2 {\±) <S> |±)} then the two-qubit graph basis measurements are 
defined as Mf 2 = AZ X2 {M\ <g> M 2 °) AZ 12 . It is now natural to extend our definition of Af-p 2 to obtain 
the family of 2-qubit measurements of interest: 

M%f := AZ 13 (Mf g> Mf) AZ 12 (17) 

corresponding to projections on the basis Q a fi '■= f\Zi 2 {P\{a) <8> P 2 ((3))({\±) (8) This family of 

two-qubit measurements together with the preparation, entanglement and corrections commands 
of the one-way model define the teleportation model. 

Before we carry on, a clarification about our choice of measurements in the teleportation model 
is necessary. The usual teleportation protocol uses Bell basis measurement defined with 

B = AX 12 {\±) ® |0/1» 
Mf 2 = AX 12 (Mf ®Mf)AXi2 

where M z is the computational-basis measurement. Note how similar these equations are to the 
equations defining the graph basis measurements. This is a clear indication that everything that 
follows can be transferred to the case where X replaces Z, and B replaces Q. However, since the 
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methodology we adopt is to embed the 2-qubit measurement based model in the one-way model, 
and the latter is based on AZ and G, we will work with the graph-basis measurements. Furthermore, 
since AZ is symmetric, whereas AX {a.k.a. as C-NOT) is not, the algebra is usually nicer to work 
with. 

Now we prove that the family of measurements in Equation [T7] leads to a universal model, 
which embeds nicely into the one-way model, but first we need to describe the important notion 
of dependent measurements. These will arise as a consequence of standardization; they were not 
considered in the existing teleportation models. 

In what follows we drop the subscripts on the AZ unless they are really necessary. We write 
(s(i),s(j)) S Z2 x Z2 to represent outcome of a 2-qubit measurement, with the specific convention 
that (0,0), (0,1), (1,0), and (1,1), correspond respectively to the cases where the state collapses 
to AZ|+ Q )|+ Q ), AZ|+ a )|- a ), AZ|- a )|+ a ), and AZ|- Q )|- a ). 

We will use two types of dependencies for measurements associated with X-action and Z-action: 

[M?>P](>>*) = M^ s <-V tf3 

(U,V) ^^'^j _ J^a+UTT,(3+VTT 

where s, t, u and v are in Z2. The two actions commute, so the equations above define unambigu- 
ously the full dependent measurement ( n >^)[M°|'^]( s '*). Here are some useful abbreviations: 

(0,0)[ M a,/3]( S ,t) . = [ M a,/3](s,t) 
(u,v) n^fa,/3] (0,0) ._ (u,v) n^a,/?j 
(0,0)[ M a,/3](0,0) . = M a,P 
M a,x ._ M a,0 
M a,y ._ M a,% 



As in the 1-qubit measurement case we obtain the following rewriting rules for the teleportation 
model: 

EijXf XfZ^Eij EX 

EijZf ZfEij EZ 

{u,v)[ M a,P^ Sj t) Z r ^ (u+r,v)^ M a,l3^s,t) MZ 
M[ M ^(s,t) Z r ^ (u,v+r){ M a,^(s,t) MZ 

to which we add also the trivial commutation rewriting which are possible between commands that 
don't overlap (meaning, acting on disjoint sets of qubits). 



8.2.1 Embedding 

We describe how to translate 2-qubit EMC patterns to 1-qubit patterns and vice versa. The 
following equation plays the central role in the translation: 

M5f = E ij (M?®MP)E ij (18) 

Note that this immediately gives the denotational semantics of two-qubit measurements as cptp- 
maps. Furthermore, all other commands in the teleportation model are the same as in the one-way 
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model, so we have right away a denotational semantics for the entire teleportation model in terms 
of cptp-maps. 

We write for the collection of patterns in the one-way model and X for the collection of 
patterns in the teleportation model. 

Theorem 6 There exist functions [•]/ : V — > T and [•]& :T^V such that 

1. VVey:lV\ = {[V\f\; 

2. VTGT:[TJ = [[T] 6 ]; 

3- [']/ ° Y\b an d ['}b ° ["]/ are both identity maps. 
Proof. We first define the forward map [•] f in stages as follows for any patterns V = (V, I, O, A n . . . A\) 

1. For any i £ V \ O (i.e. measured qubits) we add an auxiliary qubit id called a dummy qubit 
to the space V. 

2. For any i G V \ O we replace any occurrence of M" with MfM* d . 

3. We then replace each of the newly created occurrences of MfMf d by M^Ea d . 

Now we show that the first condition stated in the theorem holds; we do this stage wise. The 
first two stages are clear because we are just adding qubits that have no effect on the pattern 
because they are not entangled with any pre-existing qubit, and no other command depends on a 
measurement applied to one of the dummy qubits. Furthermore, we add qubits in the state |+) 
and measure them in the |±) basis. The invariance of the semantics under stage 3 is an immediate 
consequence of Equation [18] and the fact that all the measurements are destructive, and hence an 
entanglement command on qubits appearing after a measurement of any of those qubits can just 
be removed. 

The map is defined similarly except that there is no need to add dummy qubits. One 
only needs to replace any two-qubit measurement with MfM^Eij. Again, this clearly pre- 

serves the semantics of patterns because of Equation [18] and the above remark about destructive 
measurements. Thus condition 2 of the theorem holds. 

The fact that the two maps are mutual inverses follows easily. As all the steps in the translations 
are local we can reason locally. Looking at the forward mapping followed by the backward mapping 
we get the following sequence of transformations 

Mf ^ stag el,2 M?Mf d 

=^ Equation [18] ^ii d ^"d 
^Equation [T8l M° 'Mf d Eu d Ea d 

M?M* d 

This shows that we have the third condition of the theorem. □ 

Note that the translations are compositional since the denotational semantics is and also it fol- 
lows immediately that the teleportation model is universal and admits a standardization procedure. 
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Example. Consider the teleportation pattern in the teleportation model given by the command 



sequence: X^ 1 Z^ 2 E23, we perform the above steps: 



A 3 Z 3 



M x,x 



E- 



23 



=^ Equation Il8l 



X° 1 Z° 2 M?MZE 12 E 2 3 
and hence obtain the teleportation pattern with 1-qubit measurements. 



Example. We saw before, the following EMC 1-qubit pattern for R z (a) which can be embedded 
to an EMC 2-qubit pattern using the above steps: 



A 3 Z 3 
A 3 Z 3 



E12E23 



[M^] Sl M^M{ a M^ 



E12E23 



^ stage 1,2 

^Equation 1181 and standardization 



M^]^M- 1 a ; x E lld E22 d E 1 2E 2 3 



Note that we have explicit algorithmic translations between the models and not just illustrative 
examples. This is the main advantage of our approach in unifying these two models compared to 
the extant work [CLN051 KIMlJP05j . 



9 Conclusion 

We have presented a calculus for the one-way quantum computer. We have developed a syntax 
of patterns and, much more important, an algebra of pattern composition. We have seen that 
pattern composition allows for a structured proof of universality, which also results in parsimonious 
implementations. We develop an operational and denotational semantics for this model; in this 
simple first-order setting their equivalence is clear. 

We have developed a rewrite system for patterns which preserves the semantics. We have 
shown further that our calculus defines a polynomial-time standardization algorithm transforming 
any pattern to a standard form where entanglement is done first, then measurements, then local 
corrections. We have inferred from this procedure that the denotational semantics of any pattern is 
a cptp-map and also proved that patterns with no dependencies, or using only Pauli measurements, 
may only implement unitaries in the Clifford group. 

In addition we introduced some variations of the one-way and teleportation models and pre- 
sented compositional back-and-forth embeddings of these models into the one-way model. This 
allows one to carry forward all the theory we have developed: semantics, rewrite rules, standard- 
ization, no-dependency theorems and universality. In fact the result of making the connection 
between the one-way model and the teleportation model is to introduce ideas: dependent mea- 
surements, standard forms for patterns and a standardization procedure which had never been 
considered before for the teleportation model. This shows the generality of our formalism: we 
expect that any yet to be discovered measurement-based computation frameworks can be treated 
in the same way. 

Perhaps the most important aspect of standardization is the fact that now we can make patterns 
maximally parallel and distributed because all the entanglement operators, i.e. non-local operators, 
can be performed at the beginning of the computation. Then from the dependency structure that 
can be obtained from the standard form of a pattern the measurements can be organized to be as 
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parallel as possible. This is the essence of the difference between measurement-based computation 
and the quantum circuit model or the quantum Turing machine. 

We feel that our measurement calculus has shown the power of the formalisms developed by 
the programming languages community to analyze quantum computations. The ideas that we use: 
rewriting theory, (primitive) type theory and above all, the importance of reasoning compositionally, 
locally and modularly, are standard for the analysis of traditional programming languages. However, 
for quantum computation these ideas are in their infancy. It is not merely a question of adapting 
syntax to the quantum setting; there are fundamental new ideas that need to be confronted. What 
we have done here is to develop such a theory in a new, physically-motivated setting. 

There were prior discussions about putting patterns in a standard form [RB02] but these worked 
only with strongly deterministic patterns, furthermore one needs to know which unitary is being 
implemented. In our case the rewrite rules are entirely local and work equally well with all patterns. 

An interesting question related to the measurement calculus is whether one can give sufficient 
conditions - depending only on the entanglement structure of a pattern - that guarantee deter- 
minacy. In a related paper the first two authors have solved this problem |DK05a] . In effect 
given an entanglement structure with distinguished inputs and outputs one can enumerate all the 
unitaries that can be implemented with it. This gives a precise handle on the entanglement re- 
sources needed in the design of specific algorithms and protocols directly in the measurement-based 
model |dBDK06j . 

Finally, there is also a compelling reading of dependencies as classical communications, while 
local corrections can be thought of as local quantum operations in a multipartite scenario. From 
this point of view, standardization pushes non-local operations to the beginning of a distributed 
computation, and it seems the measurement calculus could prove useful in the analysis of dis- 
tributed quantum protocols. To push this idea further, one needs first to articulate a definition of 
a distributed version of the measurement calculus; this was done in a recent paper [DDKP05] . The 
distributed version of the calculus was then used to analyze a variety of quantum protocols and to 
examine the notion of knowledge flow in them [DP05| . 
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A Background on Quantum Mechanics and Quantum Computa- 
tion 

We give a brief summary of quantum mechanics and quantum computing. We develop some of 
the algebra, define some notations, and prove a couple of equations which we have used in the 
paper. Although the paper is self-contained, the reader will find the expository book of Nielsen and 
Chuang [NC00J useful for quantum computation or the excellent book by Peres |Per95] for general 
background on quantum mechanics. 
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A.l Linear Algebra for Quantum Mechanics 

We assume that the reader is familiar with the basic notion of a vector space. In quantum mechanics 
we always consider vector spaces over the complex numbers. For quantum computation the vector 
spaces are always finite dimensional. The vector spaces that arise in quantum mechanics are Hilbert 
spaces and are thus usually written Sj; that is they have an inner product usually written (u,v) 
where u and v are vectors. The inner product is a map from Sj x Sj to the complex numbers C. The 
inner product is linear in the second argument but aniz-linear in the first argument. In general, 
there is a topological completeness condition on Hilbert spaces but, in the finite dimensional case 
this is automatic and we will ignore it. 

Following Dirac, it is customary to call elements of Sj kets and write them in the form \u) or 
whatever symbol is appropriate inside the half-bracket. The dual vectors are called bras and are 
written (v\; the pairing thus can naturally be identified - conceptually and notationally - with the 
inner product. 

Linear operators come naturally with vector spaces; a linear operator is a linear map from 
a vector space to itself. Linear operators on finite dimensional spaces are often represented as 
matrices. The most important notion for an operator on a Hilbert space is that of an adjoint. 

Definition 17 If A : Sj — > Sj' is a linear operator then the adjoint, written A\ is a linear operator 
from Sj' to Sj such that 

Vu £ Sj', v £ Sj(u, Av) = (A^u, v). 

In terms of matrices this just amounts to transposing the matrix and complex conjugation each 
of the matrix entries; sometimes this is called the hermitian conjugate. An inner product preserving 
linear map is called a unitary embedding. When Sj = Sj' we can also define the following operators. A 
hermitian operator A is one such that A = A^ and a unitary operator U is one such that U~ l = W. 
A projection P is a linear operator such that P 2 = P and P = pt . A projection operator can be 
identified with a subspace, namely its range. The eigenvalues of a hermitian operator are always 
real. Suppose U is a unitary, and P a projection, then UPU^ is also a projection. 

It is common to use the Dirac notation to write projection operators as follows: given a vector 
\a) of unit norm, the projection onto the subspace spanned by \a) is written |a)(a|. To see why this 
makes sense, suppose that \b) is another vector then its component along \a) is the inner product 
(a, h). Now if we just juxtapose the expressions \a) (a\ and \b) we get \a)(a, b), viewing the (a, b) as 
a number and moving it to the front we get (a, b)\a) as the result, which is the right answer for the 
projection of \b) onto \a). Thus one can apply the projection operator just by juxtaposing it with 
the vector. This kind of suggestive manipulation is part of the appeal of the Dirac notation. 

One important fact - the spectral theorem for hermitian operators - states that if M is a 
hermitian operator, Aj are its eigenvalues and Pi are projection operators onto the corresponding 
eigenspaces then one can write 

M = ^2\ i P i . 

i 

If we have \i) as the normalized eigenvectors for the eigenvalues Aj then we can write this in Dirac 
notation as: 

i 

Finally we need to combine Hilbert spaces. 
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Definition 18 Given two Hilbert spaces f) with basis vectors {aj|l < i < n} and $)' with basis 
< J < w} we define the tensor product, written Sj® Sj 1 ' , as the vector space of dimension n-m 
with basis (g) 6j . 

There are more elegant, basis-independent ways of describing the tensor product but this definition 
will serve our needs. We almost never write the symbol ® between the vectors. In the Dirac 
notation this is always omitted and one writes, for example, \uv) instead of \u) (8) \v). 

The important point is that there are vectors that cannot be written as the tensor product of 
vectors. For example, we can write a\ <g> b\ + a<i <g> 62 where the ai and the i»j are basis vectors of 
two 2-dimensional Hilbert spaces. This means that given a general element of fj <g> Sj' one cannot 
produce elements of Sj and this is very different from the cartesian product of sets. This is the 
mathematical manifestation of entanglement. 

A very important function on square matrices is the trace. The usual trace - i.e. the sum of 
the diagonal entries - is basis independent and is actually equal to the sum of the eigenvalues, 
counted with multiplicity. The trace of A is written tr(A) and satisfies the cyclicity property 
tr(AB) = tr(BA); applying this repeatedly one gets 

tr{A 1 ...A n )=tr(A a[l) ...A (T{n) ) 

where a is a cyclic permutation. The explicit formula for the trace of A : V — > V is tr(A) = 
X^I-^N) where \i) is a basis for V. 

One often needs to compute a partial trace. Consider a linear map L : V (8> W — ► V <8> W. 
Suppose that \vi) is a basis for V and \ wi) is a basis for W then \viWj) is a basis for V W. Now 
we can define the partial trace over V as 

try (A) : W = ^2{vi\A\vi). 

i 

This corresponds to removing the V dependency; often we use the phrase "tracing out the V 
component." 

A. 2 Quantum Mechanics 

We state the basic facts of quantum mechanics and will not discuss the experimental basis for this 
framework. The key aspects of quantum mechanics are: 

• the states of a quantum system form a Hilbert space, 

• when two quantum systems are combined, the state space of the composite system is obtained 
as the tensor product of the state spaces of the individual systems, and 

• the evolution of a quantum system is given by a unitary operator, and 

• the effect of a measurement is indeterminate. 

The first says that one can form superpositions of the states. This is one of the most striking 
features of quantum mechanics. Thus states are not completely distinct as they are in classical 
systems. The inner product measures the extent to which states are distinct. The fact that systems 
are combined by tensor product says that there are states that of composite systems that cannot 
be decomposed into individual pieces. This is the phenomenon of entanglement or non-locality. 
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Measurement is what gives quantum mechanics its indeterminate character. The usual case, 
called projective measurements, is when the quantity being measured is described by a hermitian 
operator M. The possible outcomes are the eigenvalues of M. If M is an observable (hermitian 
operator) with eigenvalues Aj and eigenvectors \<fii) and we have a generic state \ip) = X^^I'M 
then the probabilities and expectation values of the measurement outcomes are given by: 

. Prob{\iU)) = \a\ 2 

• E[MW\ = EilctpAj = Ei M&) = (rP,M?P). 

It is important to note that the effect of the measurement is that the projection operator Pj 
is applied when the result Aj is observed. The operator M does not describe the effect of the 
measurement. 

The formulas above reveal that no aspect of a measurement is altered if the vector describing 
a quantum state is multiplied by a complex number of absolute value 1. Thus we can multiply 
a state by e %d without changing the state. This is called changing the phase. While the phase is 
not important phase differences are very important. Multiplying a vector by a complex number 
is a change of phase as well as a change in its length. Usually we normalize the state so that we 
can read the results of measurements as probabilities. Given a vector the subspace spanned by it - 
always one dimensional - is called a ray. Thus a state is really a ray rather than a vector. However, 
it is customary to blur this distinction. 

A. 3 Some qubit algebra 

Quantum computation is carried out with qubits the quantum analogues of bits. Just as a bit 
has two possible values, a qubit is a two dimensional complex Hilbert space, in other words it is 
(isomorphic to) the two dimensional complex vector space C 2 . 

One works with a preferred basis, physically this corresponds to two distinguishable states, like 
"spin up" and "spin down". One writes |0), and |1) for its canonical basis, so that any vector ip 
can be written as a|0) + f3\l) with a, (3 in C. Furthermore, C 2 can be turned into a Hilbert space 
with the following inner product: 

(a\0) +p\l),a'\0) + := a* a' + /3*(3' 

where a* is the complex conjugate of a. One then obtains the norm of a vector as: 

W\\ := <^,# = (a*a + /3*/3)i 

Given V a finite set, one writes Sjy f° r the Hilbert space <S> ue yC 2 ; the notation means an n-fold 
tensor product of the C 2 where n is the size of V. A vector in Sjy is sa id to be decomposable if 
it can be written ® u ev' t Pu for some ip u £ C 2 . Such decomposable vectors will be written e in the 
sequel. Decomposable vectors can be represented by a map from V to C 2 , and we will use both 
notations depending on which is more convenient. As we have noted before there are some vectors 
that are not decomposable. 

As in the case of C 2 , there is a canonical basis for $)y, sometimes also called the computational 
basis, containing decomposable vectors e such that for all v G V, e(v) = |0) or e(v) = |1). 

The inner product on Sjy, according to the general definition given above, is defined on decom- 
posable vectors as: 
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Note that all vectors in the computational basis are orthogonal and of norm 1. The vectors of norm 
1 are usually called unit vectors; we always assume that states are described by unit vectors as 
noted before. 

Here are some common states that arise in quantum computation: 



|o> = I T) 



ID 



I) 



,1+) 



l 

71 



i 

71 



i 



-i 



It is easy to see that a linear operator is unitary if it preserves the inner product and hence the 
norm. Thus unitaries can be viewed as maps from quantum states to quantum states. 

Some particularly useful unitaries are the Pauli operators given by the following matrices in the 
canonical basis of C 2 : 



X 



1 

1 



Y 



-% 

1 



Z 



1 

-1 



We note that all these operators are involutive, self-adjoint, and therefore unitaries. All these 
matrices have determinant = — 1 . We will not discuss the geometric significance of these operators 
here; their real importance in quantum mechanics comes from the fact that they can be used to 
describe rotations, thus they are usually called the "Pauli spin matrices" by physicists. 

Some basic algebra of these matrices are given below. First they all square to the identity. 

X 2 = Y 2 = Z 2 = I. 

The Pauli operators do not commute. We use the notation [A, B] for AB — BA, the commutator 
of A and B. The commutator measures the extent to which two operators fail to commute: it is 
customary to present the algebra of operators using it. One also uses the symbol {^4, B} to stand 
for AB + BA: it is called the anti- commutator. For the Pauli operators we have the following 
commutators and anti-commutators : 



XY = iZ YX = -iZ [X, Y) = 2iZ {X, Y} = 

ZX = %Y XZ = -iY [Z, X] = 2iY {Z, X} = 

YZ = iX ZY = -iX [Y, Z] = 2iX {Y, Z} = 

Definition 19 Define the Pauli group, P n , as the group consisting of tensor products of I, X, Y, 
and Z on n qubits, with an overall phase of ±1 or ±i. 

Given a group G the operation x i— > g~ 1 xg is called conjugation by g. These conjugations 
give the effect of switching operators around. If G is a group and H is a subgroup of G then the 
normalizer of H is another subgroup of G, say K, with the property that for all h G H , k G K we 
have k hk in H. 

The effect of conjugating measurements and other corrections by Pauli operators is a key part 
of the rewrite rules described in the main text. They can be verified using the algebra given here. 
A very important related group is called the Clifford group. 

Definition 20 The Clifford group, C n , is the group of unitary operators that leave the Pauli group 
invariant under conjugation, i.e. it is the normalizer of the Pauli group viewed as a subgroup of the 
unitary group. 
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The Clifford group on n qubits can be generated by the Hadamard transform, the controlled- X 
(CNOT) or controlled- Z (AZ), and the single-qubit phase rotation: 



H= 72{1 -1 



CNOT 



(1 0\ 

10 

1 

\0 1 0/ 



AZ 



(I \ 

10 

10 

\0 -1/ 



p 



1 

? 



The importance of the Clifford group for quantum computation is that a computation consisting 
of only Clifford operations on the computational basis followed by final Pauli measurements can be 
efficiently simulated by a classical computer, this is the Gottesman-Knill theorem [Got97l [NCOOj . 



A. 4 Density Matrices 

In order to capture partial information about quantum systems one uses density matrices. Before 
we describe density matrices we review some linear algebra in the bra-ket notation. Given a ket l^) 
the notation \ip){ip\ denotes the projection operator onto the one dimensional subspace spanned by 
To verify this note that 

and 

(h/>) <</>!) (|4>» = \t/>)(m = (</#>h/>>- 

If is an orthonormal basis for fj the identity matrix is written If Q is a linear 

operator with eigenvalues qi and eigenvectors \qi), which form an orthonormal basis for fj, we can 
represent Q as ^ Qi\Qi){li\- To see this, let \ip) = X^ c «k«) then 

% i 

now using our representation for Q we calculate 

Q\i>) = ^iMi){ii\{W)) = J2cjqi\qi}(<n\qj) = ^2,wMii- 

i i,j i 

This is a version of the spectral theorem that we mentioned in the first subsection of this 
appendix. 

A state (i.e. a ray in Sj) is called a pure state. If a and b are distinct eigenvalues of some 
observable A with corresponding eigenvectors \a) and \b) it is perfectly possible to prepare a state 
of the form 7/g(l a ) + |&))- A measurement of A on such a state will yield either a or b each with 

probability ^. However, it is also possible that a mixture is prepared. That is to say instead of 
a quantum superposition a classical stochastic mixture is prepared. In order to describe these we 
will use density matrices. 

For a system in a pure state the density matrix is just the projection operator \ifj){ip\. If 
we have an observable Q with eigenvalues qi - assumed nondegenerate for simplicity - then we can 
expand 1^} in terms of the eigenvectors by 

W) = z^cMi)- 

i 



43 



Now the probability of observing qi when measuring Q in the state \ip) is |(gi|V>)| 2 . Recalling that 
the identity is given by I = £^ ■ \qj)(qj\ we get that 

Prob(qi,\^}) = ^2(qi\4>)(4>\qj}(qj\qi) 

3 

which after rearranging and using the definition of trace of an operator yields 

Tr((\ qi )( qi \)(\m\))- 

If as is typical we write p^ for the density matrix and Pi for the projection operator onto the 
subspace spanned by the eigenvector \qi) we get 

Prob( qi , |V>)) = Tr{P iP ). 

It is an easy calculation to show that the expectation value for Q in the state \tp) is Tr(Qp). 

What if the state is not known completely? Suppose that we only know that a system is one of 
several possible states \ipi), . . . , \tpk) with probabilities pi,...,Pk respectively. We define the density 
matrix for such a state to be 

p = ^2pi\A)(^i\- 

1=1 

The same formulas for the probability of observing a value qi , i.e. Tr(Pip) and for the expectation 
value of Q, i.e. Tr(Qp) apply. One can check directly that a density matrix has the following two 
properties. 

Proposition 21 An operator p on S) is a density matrix if and only if 

• p has trace 1 and 

• p is a positive operator, which means that it has only positive eigenvalues or, equivalently, 
that for any x S Sj we have (x, px) > 0. 

Furthermore, if p is a density operator, Tr(p 2 ) < 1 with equality if and only if p is a pure state 
(i.e. a projection operator). 

Suppose that we have a density matrix p describing a pure state of an n+m dimensional system. 
Now suppose that an observer can only see the first n dimensions. The density matrix £ describing 
what he can see is contained by taking the partial trace over the m dimensions that the observer 
cannot see. Doing this gives, in general, a nonpure state. Similarly a complementary observer who 
sees only the m dimensions would construct her density matrix a by taking the appropriate partial 
trace. Taking these traces loses information; in fact, one cannot reconstruct p even from both £ 
and a. Certainly the tensor product of £ and a does not give back p. This is due to the loss of the 
cross-correlation information that was encoded in p but is not represented in either £ or a. 

The axioms of quantum mechanics are easily stated in the language of density matrices. For 
example, if evolution from time ti to time ti is described by the unitary transformation U and p 
is the density matrix for time t±, then the evolved density matrix p' for time t2 is given by the 
formula p' = UpU^ . Similarly, one can describe measurements represented by projective operators 
in terms of density matrices [NC001 IPre98j . Thus if a projector P acts on a state \tp) then the 
result is P\ip); the resulting transofrmation of density matrices is For a general 

density matrix p we have p i— > PpP, note that since P is self-adjoint we do not have to write P*. 
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A. 5 Operations on Density matrices 

What are the legitimate "physical" transformations on density matrices? Density matrices are 
positive operators and they have trace either equal to 1 if we insist on normalizing them or bounded 
by 1. These properties muct be preserved by any transformations on them. 

We need first to define what it means for a vector to be positive. Any vector space V can be 
equipped with a notion of positivity. 

Definition 22 A subset C ofVis called a cone if 

• x £ C implies that for any positive a, ax G C , 

• x,y £ C implies that x + y £ C and 

• x and —x both in C means that x = 0. 

We can define x > to mean x £ C and x > y to mean x — y £ C . 
Definition 23 An ordered vector space is just a vector space equipped with a cone. 
It is easy to check the following explicitly. 

Proposition 24 The collection of positive operators in the vector space of linear operators forms 
a cone. 

Now we can say what it means for a map to be a positive map. 
Definition 25 Abstractly, L : (V, <v) — ► (W, <w) ^ s a positive map if 

Vf £ V. v> v 0^ L(v) > w 0. 

It is important to not confuse "positive maps" and "positive operators." 

If we are transforming states (density matrices) then the legitimate transformations obviously 
take density matrices to density matrices. They have to be positive maps considered as maps 
between the appropriate ordered vector spaces. The appropriate ordered vector spaces are the 
vector spaces of linear operators on fj the Hilbert space of pure states. 

Unfortunately the tensor product of two positive maps is not positive in general. We really 
want this! If one can perform transformation T\ on density matrix p\ and transformation T2 on 
density matrix pi then it should be possible to regard p\ ® pi as a composite system and carry out 
T\ ® T2 on this system. We certainly want this if, say, T2 is the identity. But even when T2 is the 
identity this may fail; the usual example is the transposition map, see, for example [NCOOj . 

The remedy is to require the appropriate condition by fiat. 

Definition 26 A completely positive map K is a positive map such that for every identity map 
I n : C" — > C n the tensor product K ® I n is positive. 

It is not hard to show that the tensor of completely positive maps is always a completely positive 
map. This condition satisfies one of the requirements. We can insist that they preserve the bound 
on the trace to satisfy the other requirement as well. However we would like an explicit way of 
recognizing this. 

The important result in this regard is the Kraus representation theorem [Cho75j . 
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Theorem 7 (Kraus) The general form for a completely positive map £ : B(S)\) — > B{$)2) is 

m 

where the A m : S)\ — > 5)2- 

Here B(Sj) is the Banach space of bounded linear operators on Sj. If, in addition, we require that 
the trace of £{p) < 1 then the A m will satisfy 

in 

The following term is common in the quantum computation literature. 

Definition 27 A superoperator T is a linear map from 23y to *Bu that is completely positive and 
trace preserving. 



46 



